function ejer6b
        
    A = [5 7 6 5; 7 10 8 7;6 8 10 9; 5 7 9 10] 
    b = [23; 32; 33; 31]
    x = [0; 0; 0; 0];
	er = 10^-4;
	m = 10000;
      
    fprintf('solucion de jacobi \n');
    [x1,k1] = jacobi(A, b, er, m, x)
    
    fprintf('solucion de gauss seidel \n');
    [x2,k2] = gseidel(A, b, er, m, x)
    
 endfunction   